Two inequivalent sublattices and orbital ordering in MnV 2 04 by 51 V NMR 
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We report detailed 51 V NMR spectra in a single crystal of MnV2C>4. The vanadium spectrum 
reveals two peaks in the orbitally ordered state, which arise from different internal hyperfine fields 
at two different V sublattices. These internal fields evolve smoothly with externally applied field, 
and show no change in structure that would suggest a change of the orbital ordering. The result 
is consistent with the orbital ordering model recently proposed by Sarkar et al. [Phys. Rev. Lett. 
102, 216405 (2009)] in which the same orbital that is a mixture of t2 g orbitals rotates by about 45° 
alternately within and between orbital chains in the I4\/a tetragonal space group. 



Competition between spin interactions and orbital de- 
generacy is a key factor in determining the ground state 
magnetic and lattice structure of transition metal ox- 
ides, and gives rise to a rich spectrum of phase transi- 
tions in magnetic insulators. 1,2 Recently, vanadium ox- 
ide spinels of the form AV2O4 have attracted interest, 
in which A is a divalent transition element that is ei- 
ther non-magnetio 3 -^ 5 - (Mg, Zn, Cd) or magnet i(JL2&iiL12. 
(Mn) . The intriguing physics of the vanadate compounds 
arises from the geometrical spin frustration of the triply 
degenerate tig orbitals of V 3+ (3c? 2 , S = 1) in the spinel 
structure. The V ions sit at the vertices of corner-sharing 
tetrahedra, and experience a magnetic exchange interac- 
tion with nearest neighbor V spins. In the process of re- 
lieving spin frustration, the vanadates typically undergo 
two consecutive phase transitions. First, a structural dis- 
tortion splits the tig levels into a low-lying xy orbital 
and a higher doublet (yz, zx) at a temperature Tg. This 
structural phase transition is accompanied by long-range 
orbital ordering of the V orbitals, since Hund's rules im- 
ply that a single electron occupies the excited doublet 
and orbital exchange interactions lift this degeneracy. If 
the A site is non-magnetic, then the V spins order anti- 
ferromagnetically at a lower temperature TV < I5. 

When the A site is magnetic, novel features emerge 
that are distinct from the other vanadates. In partic- 
ular, for Mn 2+ (3d 5 ) a ferrimagnetic (FEM) transition 
occurs before the structural transition (Tg < TV). The 
combination of strong ferromagnetic Mn-Mn couplings, 
antiferromagnetic (AFM) Mn-V couplings, and AFM V- 
V couplings leads to a collinear FEM spin configuration 
at TV ~ 56 K. This collinear state retains the orbital 
degeneracy, but is unstable so that upon further cooling 
the V spins become non-collinear below the structural 
transition at Tg = 53 K. Once again, long range orbital 
order of the V d orbitals emerges in concert with struc- 
tural distortion. However, the nature of the orbital sym- 
metry remains poorly understood and continues to be 
debated theoretically in the literaturei n i 12 ' 13 i 14 Experi- 
mentally, Adachi et al. proposed the ferro-orbital order- 
ing in the IAi/amd tetragonal symmetry^ but later x-ray 
and neutron experiments on single crystals^i 1 ^ claimed 



that the 74 1 /a space group was present with antiferro- 
orbital ordering in which yz and zx orbitals alternate 
along the c axis. Although I&i/a space group allows only 
antifcrro-orbital order due to the symmetry considera- 
tion, Chung et ali^ 5 - pointed out that the large exchange 
coupling along c determined by neutron scattering is con- 
tradictory with the simple antiferro-orbital ordering. In 
fact, taking into account the staggered trigonal distor- 
tion, recent first principles calculations proposed a new 
orbital ordering in which the same orbital forms "anti- 
ferro" orbital order by rotating its direction by ~ 45° al- 
ternately within and between orbital chains, while main- 
taining IA\/a symmetry^ The newly proposed orbital 
ordering model reconciles the contradicting experimental 
results and, indeed, is in good agreement with our NMR 
results. The three proposed orbital ordering models are 
schematically drawn in Fig. [TJ 

In order to investigate the orbital and magnetic order 
microscopically, we have carried out 51 V nuclear mag- 
netic resonance (NMR) as a function of field (H) and 
temperature (T) within the FEM state. In zero field, the 
51 V nuclear spins levels are split by an internal hyperfine 
field from the ordered Mn and V moments. The field de- 
pendence of the resonance reveals the non-collinear na- 
ture of the ordered moments. Single ferrimagnetic do- 
main is formed above H c ~ 0.3 T, while the orbital or- 
dering is intact up to 13 T. 

51 V NMR spectra were obtained on a single crystal of 
MnV2C>4 between 4 and 35 K in zero field and in external 
fields up to 13 T. The preparation of the single crystal 
of MnV2C>4 has been described in detail in Ref. The 
spectra were obtained by integrating averaged spin echo 
signals as the frequency was swept through the resonance 
line. The V resonance in the single crystal (SC) is close 
to 285 MHz as we found in previous measurement of a 
polycrystal (PC) sample. 16 However, the spectrum in the 
SC is narrower than in the PC and reveals two sharp 
features, whereas the PC sample is significantly broader 
with poorly resolved features (Fig. [2]). These differences 
are consistent with recent x-ray diffraction and magne- 
tization measurements^ as well as specific heat data, 8 
which show evidence of an impurity cubic phase result- 
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FIG. 1: (Color online) Proposed orbital ordering models: (a) 
staggered antiferro-orbital order in 74i/a symmetry;*^ (b) 
ferro-orbital order in IAi/amd symmetry, 7 (c) "antiferro"- 
orbital order in I4\ /a symmetry^ Blue (solid) and red (dot- 
ted) lines represent orbital chains running along the edges, 
corresponding to [110], of the corner-shared V tetrahedra. 
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FIG. 2: (Color online) 51 V spectrum of single crystal at zero 
field. Compared with spectrum of polycrystal, the structure 
is much simple and narrow, indicating that single crystal is 
relatively free of the impurity cubic phase. The contribution 
of small impurity phases remaining in the single crystal are 
indicated by arrows. Two main peaks in single crystal spec- 
trum implies the existence of inequivalent V sublattices, as 
depicted in inset, in which hyperfine fields at V nuclei (empty 
arrows) are not antiparallel to the V 3+ spin moments. The 
heights of two spectra were nomalized for comparison. 



ing from non-stoichiometric crystallites in PC samples. 8 
The spectrum of SC clearly shows two narrow lines, but 
also reveals weak signals, in particular, near 270 MHz 
and 305 MHz, due to a small portion of the impurity 
cubic phase remained in SC sample. 

In the absence of quadrupolar effects, the V resonance 
frequency is given by vq = jnH^ where 77V = 11.195 
MHz/T is the nuclear gyromagnetic ratio, and ffhf is the 
hyperfine field at the V. A priori, one might expect that 
the two well resolved peaks might arise from the sec- 
ond order quadrupolar contribution in the orbitally or- 
dered state. In this case, the splitting can be written asi£ 
Au = 25^/144^ [/(/+l)-3/4] and the NQR frequency 
vq is estimated to be ~ 40 MHz. Since, however, we are 
unable to find any satellite transitions associated with 
the first order quadrupole splitting, we conclude that vq 
is negligibly small and the origin of the two peaks is not 
quadrupolar in nature. Thus the two peaks must arise 
from two different hyperfine fields at different V sublat- 
tices which exist only in a state with both long-range 
orbital and spin order. 

It is clear that there is only one V site in the cubic 
spinel structure above T3. The proposed IAi/a tetrag- 
onal structure below T$ exhibits a slight contraction of 
the VOe octahedra along the c direction, but still retains 
a single crystallographic V site However, the hyperfine 
fields at V nuclei can be differentiated. The net hyperfine 
field at the V can be written as the sum of on-site and 
transferred terms : H hf = H™" sitc + Hy™. Each term 
can be decomposed as: 



Hgr ito = H F + H ; + H d , (1) 

H*r s = E^ si +E^ s Mn> (2) 

i 3 

where Hp is the Fermi contact field arising from the 
core polarization, H/ the orbital term from the orbital 
momentum L, the dipolar field from the on-site elec- 
trons, the indices i and j are over nearest neighbors V and 
Mn spin moments, respectively. For an orbital triplet ion, 
the orbital momentum is only partially quenched— and 
could be as large as 0.34 [Ib in MnV204^ causing large 
H; and H^. On the other hand, the transferred terms 
are usually very small compared to the on-site termji 9 - 
but could be large enough to produce two peaks. 

In the case of the perfect antiferro-orbital orde r 9 ' 10 ! 11 
or ferro-orbital order,— ^ the hyperfine fields should be 
identical for all of the V sites. Therefore, our data sug- 
gest that (i) the tetragonal symmetry is lower than 74 1 / a 
to allow the two V sites, or (ii) the orbital ordering is 
somewhat complex, yet without lowering the tetragonal 
symmetry IA\/a. Case (i) is ruled out basing on the neu- 
tron and x-ray result o 9 ' 10 so that we assume that case (ii) 
is applicable. Indeed, case (ii) is realized in the orbital 
ordering model proposed by Sarkar et al»« in which the 
same orbital rotates alternately within and between or- 
bital chains by about 45°. Since the model was obtained 
by taking into consideration the trigonal distortion of 
VC>6 octahedra in addition to the tetragonal contraction, 
it naturally explains the two V sublattices still in the 
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FIG. 3: (Color online) Temperature dependence of 51 V spec- 
trum in zero field. A Boltzmann correction by multiplying T 
was made. The signal becomes weak rapidly with increasing 
T and disappears above 35 K due to the shortening of T2 as- 
sociated with thermal fluctuation. Dotted lines are Eq. ©, 
implying that the spin wave theory is applicable at low T re- 
gion. Seemingly different T dependence of relative intensities 
between two peaks is ascribed to the impurity cubic phase in 
which resonance frequencies are almost T independent.— 



I4:i /a tetragonal symmetry because the V-V and Mn- 

V transferred hyperfine fields become different for two 

V sublattices depending on the orbital overlap between 
V-V and/or V-Mn via oxygen 2p orbitals [see Fig.QJc)]. 

The two differentiated V sites appear to be robust 
against the variation of temperature and external field. 
As shown in Fig. [31 with increasing T, the spectrum shifts 
to lower frequencies, while retaining the same shape. 
Since -f/hf is proportional to the sublattice magnetization, 
we expect v(T) to follow Bloch's T 3 / 2 law at T < T N . 
If there is an energy gap, E g , in the spin wave excitation 
spectrum, v(T) at sufficiently low T is given by: 



v{T) - 1/(0) 



I - aT 3 ' 2 e- E ^ T 



(3) 



where a is a fitting parameter^ The dotted lines in 
Fig. [3J are fits to this expression using E g — 17 .4 K. 10 
The fit is excellent for T < 15 K, but deviates above 
T > 0.25T/V. Unfortunately, we lose the signal near 30 
K due to short spin-spin relaxation rates, I/T2, which 
increase as thermal fluctuations set in near the phase 
transition at TV = 56 K. 

We now turn to the external field dependence of the 
spectrum. Figure [H shows the two peaks as a function 
of H at 7 K. The initial slope of v(H) vs. H is zero 
as seen in the inset of Fig. [U With increasing H, the 
slope increases gradually and reaches a fixed value above 
H c ~ 0.3 T. This behavior indicates the existence of 
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FIG. 4: (Color online) External field dependence of resonance 
frequencies of two peaks at 7 K. Low field data are enlarged in 
inset. Zero slope at H = indicates the existence of domain 
structure. Excellent fit above 0.3 T to 13 T with Eq. (0} was 
obtained with the fixed set of parameters (see text) , indicating 
the formation of single domain above 0.3 T, and the robust 
spin structure up to 13 T. Inset shows the magnetization M 
in field along [110] direction at the same T. M is saturated 
to 3.2 fiB as expected from the non-collinear spin structure. 



the FEM domain structure in zero field and its align- 
ment along the external field H , forming a single domain 
above H c . This is also consistent with the magnetization 
M in field that is saturated to 3.2 he at H c , as shown 
in the inset of Fig. [U Interestingly, M becomes zero as 
H — > 0. The absence of the remanent field, despite the 
"hysteresis" in field, can be understood by assuming the 
strong coupling between the structural and the magnetic 
domains. In field, the two domain structures can be de- 
coupled because the directions of the tetragonal domains 
are limited to three crystallographic axes while the mag- 
netic ones can be aligned in any direction. Therefore, our 
data suggest that the aligned tetragonal domains^ along 
H become random as H — > 0, and so do the magnetic 
ones due to the strong coupling between them. 

The resonance frequency of the V is given by the net 
vector sum of the hyperfine and external fields at the 
nucleus: 



u(H) = 7JV |H hf + H| 



In 



H 2 



2H hf Hcos9, 
(4) 

where 9 is the angle between H^f and H. This equation 
is plotted in Fig. [4] as a solid line, and clearly fits both 
peaks up to 13 T with the two parameters Hm and 9 for 
each peak. We find that the fields and angles are 24.94 
T with 9 = 27.5° and 25.23 T with 9 = 28.8°. Since 
Mn 2+ is an orbital singlet ion (3d 5 , L = 0), the Mn 2+ 
moment can be easily aligned along H due to its large 
moment^ (4.2 and the essentially isotropic suscep- 
tibility. Then the Mn-V and V-V exchange interactions 
maintain the non-collinear FEM spin structure, which is 
depicted in the inset of Fig. [21 The excellent fits of data 



4 



to Eq. (|4]) indicate that the non-collinear spin structure 
remains robust up to 13 T. We emphasize that the dif- 
ferent angles and hyperfine fields at two V sites are in- 
consistent with either simple antiferro- or ferro-orbital 
ordering because the V site are equivalent in both cases. 
Again, however, "antiferro" orbital chains formed by ro- 
tating orbitals with respect to each other— is compatible 
with the two inequivalent V sites because the hyperfine 
fields at 51 V nuclei can vary depending on the direction 
of the orbitals and/or on different overlaps of the on-site 
orbital with surrounding V d or O 2p orbitals. 

We note, however, that the angles we obtain (~ 28°) 
differ significantly from those measured by neutron scat- 
tering (~ 65°)£ii£ This result suggests that either (i) the 
orientation of the ordered spins changes in field, or (ii) 
Hhf is not coincident with the direction of the ordered 
V moments. Case (i) is possible if the delicate balance 
of exchange fields that gives rise to the particular FEM 
structure in zero field is modified by the presence of an ex- 
ternal field. However, case (i) seems unlikely because the 
magnetization data at 7 K shows the saturation moment 
of about 3.2 hb (inset of Fig. 2|), which is very close to 
the expected value of 3.1 [Ib from the given moments and 
the angle of ~ 65°. 10 Therefore, the large tilted angles of 
~ 35° between — S and Hhf requires that the on-site or- 
bital and dipolar terms be the same order as the isotropic 
Fermi contact term that is parallel to — S. For the V 3+ 
ion, we can estimate the Fermi term Hp ~ —28 Ti21 
The anisotropic dipolar term can be estimated from the 
relation H d = 4/7(r- 3 )/! B = 2/7 x 125(r- 3 ) a . u . kG^ 



Using (r" 3 ) a . u . = 3.217^ H d ~ 11.4 T. For the or- 
bital term, the magnitude of Hi could be approximated 
as 125(r _3 ) a . u , = 40 T for fully unqucnchcd angular 
momentum^ Taking into account the quenching, Hi is 
expected to be the same order as the dipolar term, and 
is not necessary to be parallel to Hp, resulting in the 
total hyperfine field that is quite off the direction of the 
ordered moment. In this sense, the large tilted angle of 
Hht could be understood as a consequence of the com- 
plex orbital ordering in the orbital triplet ground state 
of V 3+ ion, since the angular momentum is quenched in 
the perfect antiferro orbital ordering. 9 

In conclusion, our 51 V NMR study on a single crys- 
tal of MnV2C>4 reveals the two inequivalent V sublat- 
tices that imply a complicated orbital ordering pattern 
in the IA\/a symmetry. Although we cannot determine 
a specific orbital ordering solely by NMR, our data put 
a strong constraint on the possible orbital ordering mod- 
els. We find that the model proposed by Sarkar et alJ^ 
is quite promising being compatible with our data. Also 
we have shown that the orbital ordering and the non- 
collinear spin structure are robust up to 13 T. 
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